Numerical modelling approach for estimation of a yield zone in the face of a deep longwall panel

In the Longwall mining method, the coal face is being supported with two-legged or four-legged shield supports which improve better roof control between the face and canopy tip. The development of load on the face increases with the face retreating rate or the increase of the overhung length behind the shield support and as a result, the coal face tends to yield/fail. The yield/failure zone in the face extends depending on intensity of developed load due to the depth of mining, panel size, height of excavation and geometry of the overhung. The broken coal from the yielded face may fall on the shearer, armoured flexible conveyor (AFC), shield supports or may hit the deployed manpower in the face and causes the face stoppage which in turn results in to the loss of the project. It is therefore, a study is conducted to understand the face yield/failure behaviour considering various geo-mining conditions based on input parameters such as main roof thickness, overhung length, and material type and setting load. This paper develops a unique statistical model to predict the yield/failure zone in the longwall face. For this purpose, a total of fifty-four (54) three-dimensional finite element models in ANSYS software are developed and analyzed considering Drucker-Prager failure criterion.

Bai et al. 16 investigated the coal wall spall in the longwall face using FLAC2D software and pointed out that the maximum vertical stress of 38.5 MPa is developed at 3.2 m in front of the face and then it declines very sharply at 1 m of the face.Also, the horizontal deformation of the wall was observed as a maximum when the face was at less than 2 m from the monitoring station.Aghababei et al. 4 examined the longwall face floor failure using risk analysis for Parvedh-I coal mine.Prusek et al. 17 investigated the roof fall risk in longwall mines using a novel algorithm combining the empirical methods and expert opinion and reported that a yield zone in the face and cavity in the roof of 2 m and 3 m respectively.Song et al. 18 investigated the stability of longwall face using 2D finite element modelling technique and reported that maximum of 2 m yield zone occurs in the immediate roof and floor.Song and Chugh 19 analyzed the stability of longwall face deployed in the thick seams and observed that the maximum yielding of 2 m in the face and also mentioned that the structural failure ahead of the longwall face depends on the seam height, properties of coal, immediate roof and floor, loading characteristics of gob and performance of powered supports.Kong et al. 20 investigated the various parameters such as mining depth, mining height, cohesion of rock mass and support strength effecting the failure at the longwall face using 2D Universal Distinct Element Code software.Among these parameters, mining depth, mining height and cohesion showed significant influence on the face failure.Tien et al. 21conducted the study to determine deformation and failure of the coal wall at longwall face by developing physical model.Wojtecki et al. 22 carried out the study for the influence zone of the failure in overlying and underlying rock from the seam under extraction using seismic hazards.Murmu and Budi 23 investigated the longwall face failure during the main weighting period by estimating six parameters such as front abutment pressure, horizontal displacement, yielding condition, cohesion, plastic shear strain and plastic tensile strain, and reported that maximum yield zone of 5 m occurred.
From the above, it is understood that few studies were undertaken for estimation, analysis, and prediction of the yield zone of the face.Hence, it is imperative to develop a procedure to determine and predict the yield zone of the longwall face to extract the coal safely.Therefore, in this study, fifty-four (27 for setting load condition and 27 for web cutting condition) three-dimensional numerical models were prepared considering different variations in main roof overhung such as length (MRL), thickness (MRT) and material type (MRM).The yield zone data of all the numerical models were extracted and analysed to develop a novel index namely, yield zone index (YZI) to forecast the yield zone of the longwall face.

Longwall mine site description
Adriyala Longwall Mine is located in the ALP area of the Singareni Collieries Company Limited, Telangana, India.The mine has three seams such as 1, 2 and 3 having an average thickness of 6.74 m, 6.1 m and 10 m respectively.The longwall technology is deployed in the seam 1 and two panels (1 and 2) were excavated successfully and panel 3 is under extraction.This paper is prepared based on the study conducted for panel 1.The depth of the study panel and mining height were 417 m and 3.2 m respectively.The shield supports of 1.71 m wide were deployed in the longwall face of 250 m width.The length of the panel was 2500 m.The maximum capacity of the shield support was 2 × 1152 T.
Figure 1 shows the borehole section of the panel 1.From the thickness of seam 1 (6.74 m), bottom 3.2 m was excavated by the longwall and the rest thickness of 3.54 m was coal, clay and thin streaks of shaley coal considered as immediate roof.The roof rock strata lying above the immediate roof was sandstone (variable in thickness along panel length) considered as main roof and the rock layer below the coal seam/floor of the working was also sandstone.Therefore, the shield supports interact with three types of rock strata such as coal, clay and sandstone.
The first main fall in the panel 1of the Adriyala mine was reported at a retreat distance of 25 m from the setup room and the subsequent periodic falls were noticed between 10 and 15 m.It was also observed that the Physico-mechanical properties of the Adriyala mine site such as compressive strength (σ ci ) , modulus of elasticity (E ci ) and geological strength index (GSI) were collected and given in Table 1.It was noticed that the compressive strength and elastic modulus of sandstone (main roof) vary along the panel length as well as its thickness.As given in Table 1, m i is Hoek-Brown parameter and D is disturbance factor was considered as 0 as the longwall method does not have any blasting effect 2 .
The intact rock properties as mentioned in Table 1 were processed in RocLab software to estimate the rock mass properties of study panel.Table 2 shows the rock mass properties such as compressive strength (σ cm ) , modulus of elasticity (E cm ) , cohesion ( C m ), friction angle ( φ m ), dilation angle ( δ m ) based on 24 , density (ρ) and poison's ratio (ν).
After detailed investigation of the boreholes in the study panel, it was found that thickness of main roof overhung (MR) lies between 10 m and 22.5 m at an average of 15 m.Hence, in the study, MRT was considered as 10 m, 15 m and 22.5 m respectively as shown in Table 3.As mentioned earlier, the periodical falls were observed between 10 m and 20 m at an average of 13 m.Therefore, MRL was taken as 10 m, 13 m and 20 m.
Table 3 lists the variations of parameter of main roof overhung (MR) and Table 4 shows the MR overhung material type of Table 3.

Numerical model of longwall panel consists of shield support and other rock strata
A full-scale three-dimensional numerical model was prepared based on the borehole section of the study area (as mentioned in Fig. 1), panel data and engineering diagrams of the shield support.Each model consists of shield support, coal seam, immediate roof, main roof, coal bearing strata and goaf.The shield support of 2-legged 1152 T capacity was modelled and shield contains canopy, goaf shield, lemniscates links, base and hydraulic legs.The size of the numerical model was 1.71 m (equivalent to the width of the shield support) wide, 390 m length and 99.66 m height.The height and depth of mining were 3.2 m and 417 m respectively.Figures 2 and 3 show the development of 3D numerical model of the longwall panel.The clearance between the face and canopy tip was 0.1 m.The thickness of rock layer below the floor and above the roof of working was considered as 25 m and 71.46 m respectively.The additional pressure of 7.276 MPa or 345.48 m overburden was applied on top of the model.As shown in Fig. 2, the faces perpendicular to X-axis and Z-axis were applied horizontal stress (σ h = 1.5 × σ h ) and constrained in Z directions respectively.Bottom of the model was also constrained in Y-direction.
Fifty-four (54), 3D numerical models were developed considering the main roof overhung parameters such as MRL, MRT and MRM.Each parameter has three variations and these variations were done for both setting load and web cutting conditions.For simulating setting load models, the equivalent force of setting pressure (P s ) was 0.6 times of yield pressure (45 MPa) or 27 MPa was applied and estimated the setting strain.This strain was used for the development of multi-linear stress strain relationship and applied for hydraulic legs in web cutting models.All the condition of setting load models was kept same in web cut models, however, a web cut of 0.85 m was deleted and hydraulic legs in the form of bar element were inserted.All the numerical models were analysed based on the Drucker-Prager yield criterion.Detailed setting and wet cutting conditions were given in 25 .

Vertical stress distributions
Vertical stress intensity factor was defined as the ratio between induced stress and insitu stress during the retreating of face and its profile was estimated just above the shield support/roof along the panel length.Figure 4 shows the stress intensity factor for setting load (dotted line) and web-cutting conditions (dark line).From this Figure, it was observed that high tensile stress develops at the tip of the canopy and compressive stress develops on the rear side of the canopy.As expected, that high stress developed in the web cutting models due to cutting of 0.85 m web from the face.

Roof convergence
The roof convergence in the longwall face was obtained after subtracting the vertical displacement of the setting models from the that of web cutting models 26    www.nature.com/scientificreports/

Deterioration of the longwall face area
Vertical displacement distributions around the shield and longwall face area are shown in Fig. 6.It is clear that the coal face is yielded and low stress is developed near the face area.This implies that dilation of coal face is occurring and is evident from the vertical as well as the horizontal displacements as shown in Fig. 7.This Figure shows that the coal face bulges towards the shield.This phenomenon was observed in all the models having different amount of dilation depends on main roof overhung length (MRL), overhung thickness (MRT) and overhung material type (MRM).

Methodology for estimation of the yield zone in the longwall face
As mentioned above, clearance between the tip and the longwall face was kept as 10 cm in all the finite element models.After cutting of the web cut of 0.85 m, the tip to face clearance became 0.95 m.For estimation of the yield zone in the longwall face, the maximum plastic intensity of 0.01 was kept for all the numerical models.
Then, the length of the yield zone from the face line to inside the longwall face was measured as ' l1 ' for setting load models as shown in Fig. 8a.Similarly, for web cutting models, the yield zone length was measured from the face line as ' l2 ' (Fig. 8b).The yield zone observed in the setting load models ( l1 ) was found as almost the same as the yield zone of the web cutting models ( l2 ).Hence, the yield zone obtained in setting load models was not considered in further analysis.

Development of yield zone in the longwall panel
In the study, longwall face, the layers of immediate roof and a part of main roof (for a few cases) were observed as yielded.The study finds that the maximum failure zone occurs in the coal seam and extends inside the face.
From the investigation, it was found that the face and immediate roof yields severely in in all the models.The minimum and maximum yield zones were observed as 1.5 m and 3.0 m respectively.This similar phenomenon was also observed in the mine and discussed in detailed in the below section.Yield zone was estimated based on the equivalent effective plastic strain developed in the coal face and other part of the roof strata.It may be noted that higher value of effective plastic strain indicates the severe yielding.In Fig. 9a-c  It is clear that if main roof thickness is 22.5 m or T3 models, it yields irrespective of variation in other parameters.Another significant observation was that immediate roof mostly yields in all models.The part of coal seam layer exists above the main roof was also yielded in all models as shown in the Fig. 9.

Effect of the main roof overhung length
It is observed that the load on the face and shield support increases with the increment of the main roof overhung, as a result, the failure in the longwall face can be noticed and yield zone increases with the increment of main roof overhung length.It is found that growth of 11% yield zone occurred with the increment of MRL from 10 to 20 m.The yield zone of 2.51 m, 2.65 m and 2.83 m was observed for the main roof overhung length of 10 m, 13 m and 20 m respectively as shown in Fig. 10a.

Effect of the main roof overhung thickness
As mentioned in many literatures, the formation of thicker roof behind the shield support/longwall face enhances the load on the face and causes the increment of yield zone in the longwall face.In this study also, the same phenomenon was observed (Fig. 10b).For a main roof overhung thickness of 10 m, 15 m and 22.5 m, the yield zone of 2.83 m, 2.90 m and 3.0 m was observed respectively.

Effect of the main roof material type
The stiffer/strong material causes a low deformation and develops low-stress concentration on the workings than the weaker/softer rock strata.As a result, the low failure/yield zone is observed in the longwall face.From the Fig. 10c, it is also proven that the more yield zone occurred in the M1 (or softer) models, and low yield zone occurred in M3 (or strong) models.The yield zone of 2.85 m, 2.60 m and 2.50 m occurred for M1(softer), M2 (medium) and M3 (strong) models respectively.
It can be clear that in all the combinations, the immediate roof develops yield zone and may tend the cavity formation in case of L3T3 combinations of the models.This may result into falling of the roof rock pieces and face slabbing.

Validation of the yield zone observed in the numerical analysis
In this section the yield zone observed in the longwall face from all the numerical models is validated with the field observed values and the earlier published literature.
As mentioned in "Longwall mine site description" section, the average periodic weighting interval was observed between 12 and 20 m.After occurrence of periodic fall, the load on the shield support and the face was observed minimum due to occurrence of the main fall.Few field measurements of the yield zone were taken in the face with an extensometer at 4.5 m, 9.0 m, 13.5 m and 18 m face advance after a periodic fall.The yield zone of 1.62 m, 1.85 m, 2.12 m and 2.35 m were observed at face retreat of 4.5 m, 9.0 m, 13.5 m and 18 m respectively.Figure 11a-d   for cleaning the coal which resulted for delay in production, increase of the roof weighting and a reduction of periodic weighting interval.
This measured yield zone in the face is almost correlating with that of the numerical modelling results.
Very few studies [17][18][19][20]23 were undertaken to estimate the yield zone in the longwall face, floor and immediate roof. Fiure 12 shows the validation of the yield zone observed in this study with the published literature.From this Figure, it can be understood that the maximum of 5 m yield zone occurred in 23 and that of 2 m minimum occurred in 18 .It is also reported that the 79.4% of face falls extended less than 1 m, 20% lies between 1 and 2 m, 0.6% lies beyond 2 m at longwall face (Face # 1306, at a depth of 637 m) of Zhaozhuang coal mine in China.Kong et al. 20 reported that a maximum of 2.1 m face failure occurs and also mentioned that face failure extension varies proportional to cohesion of coal and support strength and inversely proportion to mining depth and mining height.It is found that an average of 2 m yield zone occurred in all other studies.However, in the present study develops the maximum yield zone of 3 m which is almost matches with them.

Development of yield zone index for prediction of the yield zone in the longwall face
As shown in the Table 5, a total of 54, 3D finite element models were developed to estimate the occurrence of yield zone in the longwall face by varying the main roof parameters such as its length, thickness and material properties.As shown in the Fig. 10a-c, the occurrence of yield zone can be noticed as an increasing trend with the growth of longer and thicker main roof and decreasing trend with the increment of material properties (from softer to strong).The variations of the yield zone with the variable parameters observed in a linear trend.As mentioned above, the yield zone in the setting load model and web cutting model was found almost same.Hence, the yield zone values (27 observations) of setting load models were not considered.Table 5 shows the yield zone values of web cutting models (27 observations) and considered in the analysis.Also, for the development of the statistical model for yield zone index, the elastic modulus of the main roof was taken to consider the effect of the material type.Figure 13 shows the 3D surface plot of yield zone occurred in the longwall face for M1 material type.www.nature.com/scientificreports/Hence, in order to combine the effect of all the input parameters (MRL, MRT and MRM), Yield Zone Index (YZI) was developed using statistical analysis to predict the yield zone for similar longwall mining condition.
In this study, YZI is the dependent variable, and input parameters such as main roof overhung length, thickness and material type are independent variables.From the statistical analysis, the best fit statistical model was developed to determine the condition of the longwall face by estimating yield zone in the longwall panel.The coefficient of determination for the model was 0.77.The yield zone index of the longwall panel was developed as given in Eq. (1).
The coefficients of the above statistical model were estimated by performing the multiple regression analysis as p 0 , p 1 , p 2 and p 3 are evaluated to be 2.321, 0.051, 0.017 and -0.112 respectively.Hence, equation of YZI is given below.
where,YZI is a yield zone index, MRL is the main roof length (m), MRT is main roof thickness (m), and MRM is main roof elastic modulus.

Conclusion
This paper proposes the methodology to estimate the yield zone of longwall face using 3D numerical models and develops the novel index "yield zone index (YZI)" by combining the effect of the main roof overhung (length, thickness and material type) to predict the condition of the longwall face before taking up the actual operation in similar geo-mining conditions.
From this study, it was observed that the yielding of the longwall face reduces as the main roof overhung becomes more competent from soft (MRM1) to strong (MRM3) and increases with increase of main roof overhung length from MRL1 to MRL3 and main roof thickness MRT1 to MRT3.The setting pressure did not show any effect on the phenomena of yielding of longwall face.
This study reveals that the maximum yield zone of 3 m occurs in the L3T3M1 model and the minimum yield zone of 1.5 m occurs in the L1T1M3 model.The results of the numerical models were also in agreement with field measured data.

Figure 1 .
Figure 1.Borehole section of the panel.

Figure 2 .
Figure 2. 3D numerical model of longwall panel consisting of shield support, goaf and other coal bearing strata.

Figure 3 .Figure 5
Figure 3. Zoomed view of 3D numerical model showing the shield, immediate roof and main roof.
. The roof convergence profiles are plotted for MRT2M2 (L1-L3) models by varying the main roof overhung lengths MRL1 = 10 m, MRL2 = 13 m and MRL3 = 20 m and they are shown in Fig. 6.A maximum vertical displacement of 42.72 mm, 43.62 mm and 48.7 mm occurred due to web cutting for main roof overhung lengths MRL1, MRL2 and MRL3 respectively.It is clearly shown that the vertical displacement on top of the coal seam near the longwall face increases with an increase of main roof overhung length from MRL1 or 10 m to MRL3 or 20 m.

Figure 4 .
Figure 4. Vertical stress intensity factor profile in the longwall panel.

Figure 5 .
Figure 5. Distribution of load on the canopy of MRL1M1(T1 to T3) models.

Figure 6 .
Figure 6.Roof convergence distribution for variation of main roof overhung length.

Figure 8 .
Figure 8. Determination of yield zone in setting load and web cutting models.
shows the yield zones of the longwall face at advance of 4.5 m, 9.0 m, 13.5 m and 18 m respectively.From these Figures, it can also be observed that most of the coal is slabbing, roof is yielding and boulders are fallen from the face.Due to this phenomenon, many time the AFC and shearer were engaged

Figure 9 .
Figure 9. Development of yield zone for various combination of main roof.

Figure 10 .
Figure 10.Behaviour of yield zone with different main roof parameters.

Figure 11 .
Figure 11.Yield zone of the longwall face observed in the study mine.

Figure 12 .
Figure 12.Validation of the yield zone observed in this study literature.
Rock strataσ ci , MPa GSI E ci , GPa D m i

Table 2 .
Rock mass material properties used in the study.

Table 3 .
Variations of the main roof (MR).